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Abstract 

We consider the gluino-mediated SUSY contributions to B° — B° mixing, B — > J/i[)K s and 
B — > Xdj in the mass insertion approximation. We find the (LL) mixing parameter can be as 
large as \(Sf 3 )LL\ < 2 x 10 _1 , but the (LR) mixing is strongly constrained by the B — > 
branching ratio and we find \{5\^)lr\ ^$ 10~ 2 . The implications for the direct CP asymmetry in 
B — > Xrf7 and the dilepton charge asymmetry (An) are also discussed, where substantial deviations 
from the standard model (SM) predictions are possible. 
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I. 



INTRODUCTION 



Recent observations of large CP violation in B — > J / ipK s |T|, |2| giving 



sin 2f3 = (0.79 ±0.10) 



(1) 



confirm the SM prediction and begin to put a strong constraint on new physics contributions 
to B° — B° mixing and B — ► J/ipK s , when combined with Ams^ = (0.472 ± 0.017) ps _1 [|J. 
Since the decay B — ► J/ipK s is dominated by the tree level SM process 6 — > ccs, we expect 
the new physics contribution may affect significantly only the B° — B° mixing and not the 
decay B — > J/ipK s . A model independent study of B° — B° mixing has been discussed 
recently by Laplace et al. [|J. However, in the presence of new physics contributions to 
B° — B° mixing, the same new physics would generically affect the B — > process. 
The relation between the new physics effects on the B° — B° mixing and B — > Xdj is in 
principle independent of each other and one may adopt a model independent analysis based 
on effective lagrangian with dimension 5 or 6 operators (for example, see Ref. || for the 
model independent study of the correlation between B — ► X s 7 and B — > X s l + l~ . The 
second paper in Ref. f| includes a new CP violating phase in the CV 7 Wilson coefficient.). 
This approach would introduce 4 new independent complex parameters: two in the B° — B° 
mixing, and two in the B — > X<ff. Being with too many independent parameters, one would 
not be able to get definite predictions based on this approach. 

In this work, we do not attempt a completely model independent study with too many 
new independent parameters. Instead, we consider B° — B° mixing, B — > J j ipK s and B^ — > 
X,ff in general SUSY models where flavor and CP violation due to the gluino mediation 
can be important. The chargino-stop and the charged Higgs-top loop contributions are 
parametrically suppressed relative to the gluino contributions, and thus ignored following 
Ref. 0. (See however Refs. 0, |j for including such effects. Another popular approach which 
is orthogonal to our approach is the minimal flavor violation model, which is discussed in 
Refs. H in the context of B physics.) We use the mass insertion approximation (MIA) for 
this purpose. Comprehensive work has been done for the first two observables in the MIA 
considering Am^ and sin 2(3 constraints only (see Ref. for the most recent studies with 
such an approach). In our work, we also include the dilepton charge asymmetry An and the 
Bd — > Xfff branching ratio constraint extracted from the recent experimental upper limit 
on the B — > fry branching ratio 



and rederive the upper limits on the (5f 3 )LL and (Sf 3 )LR mixing parameters assuming that 
only one of these gives a dominant SUSY contribution in addition to the standard model 
(SM) contribution. In addition we study the direct CP asymmetry in B^ — > X^ on the 
basis of our result for the SUSY contribution, and discuss how much deviations from the 
SM predictions are expected. Although we confine ourselves here to the gluino-mediated 
SUSY constributions only, our strategy can be extended to any new physics scenario with a 
substantial constribution to B° — B° mixing and B —> X^. 



B(B -> pj) < 2.3 x 10 
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II. EFFECTIVE HAMILTONIANS FOR B° - B° MIXING AND B X d j 



A. Effective Hamiltonian for B — B° mixing 

The most general effective Hamiltonian for B° — B° mixing (AB 
the following form ||: 



2) can be written in 



(2) 



where the operators Q^s are defined as 

Qi = 

Q 2 = 

Q 3 = 

Q 4 = 

Q 5 = 



d a R b a L <g&g 
d a A 4b% 



(3) 



and the operators Qi are obtained from Q^s by the exchange of L R. a,/3 are color 
indices, and = (1 =1=75)5/2. The Wilson coefficients Cj's receive contributions from both 
the SM and the SUSY loops: d = Cf M + Cf VSY . 

In the SM, the t — W box diagram generates only contribution to the operator Qi, and 
the corresponding Wilson coefficient Cf M at the m t scale is given by JTT| 



cfV*) 



^M* w (V; d V tb ) 2 S (x t ), 



where 



S (x t ) 



Ax t — llXt + x\ 3x| In x t 



4(l-x t ) ; 



2(1 -x t ) 



(4) 



(5) 



with Xt = m\/m\r. Performing the RG evolution down to m b scale incorporating the NLO 
QCD corrections |I2] , we get Cf M at m b : 



G 



Cf M (m b ) = ^Ml{V: d Vt b )% B S {xt)[a s {m b )]-^ 



a s (m b ) 

1 + ^r J " 



(6) 



The definition of J5 can be found in Ref. and we use the value of the QCD correction 
factor t]2B = 0.551 therein. Evaluating the matrix element of Qi, we set the bag parameter 



Bi(m b ) in the MS(NDR) scheme to 0.87 JlJ], which is numerically equal to the value in the 
RI-MOM scheme in Eqs. flTJ). 

If the deviation of the squark mass matrix from universality is small, the SUSY contribu- 
tion from the gluino-squark loop is best studied in the mass insertion approximation, which 
renders the flavor structures of the processes manifest. Flavor violations in the squark sector 
are parameterized by sizes of the off-diagonal mass matrix elements relative to the average 
squared squark mass, 



(m d ij)AB/m 



(7) 
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where i and j are family indices and A and B are chiralities, L or R. The mass matrix 
is understood to be in the super-CKM basis so that the quark-squark-gluino interaction 
vertex preserves flavor. This method is applicable to a model-independent study of flavor 
and/or CP violation in the squark sector when the series expansion in terms of (<^)ab is 
meaningful. In the presence of general (but small) flavor mixings in the down- type squark 
mass matrix, the squark-gluino box diagrams give the Wilson coefficients pi, 



SUSY 



216m 2 



24x/ 6 (x)+ 66/ 6 (x)) (S d 13 ) 



' LL 



£<SUSY 



£<SUSY 



£<SUSY 



" : 2QAxf 6 (x) (5 d 13 ) 2 



216m 2 

2 



RL 



a 



216m 2 

.,2 



36xf e (x) (S d l3 ) 



2 

RL 



504x/ 6 (x)- 72/ 6 (x)) (S d 13 ) LL (S d 13 ) 



a 



SUSY 



216m 2 

132/ 6 (x) {6 d 3 ) LR (5f 3 ) 

'2Axf 6 (x) + 120f 6 (x] 



'RR 



' RL 



at 



216m 2 



RR 



The other Wilson coefficients Cf=iJ 3 s are obtained from Cf=iJ 3 by exchange of L 
The loop functions f%(x) and fe(x), evaluated in terms of x = m^/rh 2 , are given by 



(8) 
R. 



/e(^) 



6(1 + 3x) lnx + x 3 - 9x 2 - 9x + 17 
6(x- l) 5 ' 

6x(l + x) lnx — x 3 — 9x 2 + 9x + 1 
3(x- l) 5 ' 



(9) 



These Wilson coefficients are calculated at ~ ~ m, and evolved down to the m& scale. 
A complete NLO RG evolution formula of these Wilson coefficients is available in Ref. [0. 
The initial condition (|8|) is at LO in a s , but it would be no problem to include the NLO 
correction. For this we use 



C r {mT le ) = E E {t s) + V c^) V ai C S (M S ), 



(10) 



where the SUSY scale is defined by Ms — (rh + m§)/2, and r\ = a s (Ms)/a s (mt). The list of 
'magic numbers' a.j, b[ r ' s \ and of , in the RI-MOM scheme, can be found in Ref. @. RG 
running of Ci- 3 is done in the same way as for C\- 3 . 

Each matrix element of the AB = 2 operators in (Q) is taken to be a product of its value 
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in vacuum insertion approximation and the corresponding bag parameter: 



(B d \QM 

(B d \QM 
{B d \QM 



2 



|S°) 
IS 5 ) 

IS°> 
is 5 ) 

IS 5 ) 



3 m B J Bi Bx(ji), 



5 

12 



m B„ 



1 

12 
1 
2 
1 
6 



m B d 

m b (ii) + m d (ii) 

m B d 

m b (fx) + m d {ji) 
m B . 



m 



bJbBM, 

UlBM- 



(11) 



Here we use the lattice improved calculations for the bag parameters in the RI-MOM scheme 
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SxK) = 0.87(4)+*, B 2 (m b ) = 0.82(3)(4), 
S 3 (m 6 ) = 1.02(6)(9), S 4 (m 6 ) = 1.16(3)+ 5 7 , 



\+22 
-7 ' 



(12) 



S 5 (m 6 ) = 1.91(4): 

In addition, we use the following running quark masses in the RI-MOM scheme : 

m b (m b ) = 4.6 GeV, m d (m b ) = 5.4 MeV. (13) 

The bottom quark mass is obtained from the MS mass m^' (m^ ) = 4.23 GeV. For the 
Bd meson decay constant, we assume f Bd = 200 ± 30 MeV. 

The above AS = 2 effective Hamiltonian will contribute to Am B , the dilepton charge 
asymmetry and the time dependent CP asymmetry in the decay S — > J/ ipK s via the phase 
of the S° — S° mixing. Defining the mixing matrix element by 



M 12 (S°) = 



1 



2m B 



(B°\H* B =*m 



(14) 



one has Am# d = 2|M 12 (S°)|. This quantity is dominated by the short distance contributions, 
unlike the Am^ for which long distance contributions would be significant. Therefore the 
data on Am^ p will constrain the modulus of Mi 2 (S^). On the other hand, the phase of the 

S° — S° mixing amplitude M 12 (S°) = exp(2i/3 ) |M 12 (S°)| appears in the time dependent 
asymmetry : 



A^{B° -> J/ipK s ) = sin 2/5 sinAm Sd t. 



(15) 



Since there may be large new physics (SUSY in this work) contributions to both K° — K° 
and S° - S° mixings, the CKM fit may change accordingly. Only those constraints that 
come from semileptonic processes may be used, since these will be dominated by the SM 
contributions at tree level (unless one considers R-parity violation). Therefore the angle (5 
need not be the same as the SM angle (3{— <pi), and the angle 7(= fa) should be considered 
as a free parameter in the full range from to 2ir in principle. This strategy was also adopted 
in some earlier work 0, |15|j . 
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TABLE I: Input values for the parameters. 



m Bd 
m t 

\Vcb\ 

\v ub \ 

a a (M z ) 



5.279 GeV 
174 GeV 
(40.7 ± 1.9) x 10~ 3 
(3.61 ± 0.46) x 10~ 3 
200 ± 30 MeV 
0.119 



Finally, the dilepton charge asymmetry An is also determined by Mi 2 (.B ), albeit a pos- 
sible long distance contribution to Tsm(B°). Defining the mass eigenstates of the neutral 
B° mesons as 

1 



i + M 



\B°) ±T]\B°) 



with rj = a/ (M* 2 — iY\ 2 ) / {M\2 — ^12), we can derive the following relation: 



Ai 



N{BB) -N(BB) 


V 


4 - 1 _ Im(r 12 /M 12 ) 


N(BB) + N(BB) 


V 


4 + l 1 + 


rVM 12 


2/4 



Im(r 12 /M 12 ) 



(16) 



Here Mi 2 , Ti 2 are the matrix elements of the Hamiltonian in the (B°, B°) basis: 

(B\H m \B) = M 12 - l -T l2 . 
2m B 2 

In the SM, the phases of M 12 and T 12 are approximately equal and 



AM 



SM 



2\M: 



SM 
12 h 



Ar 



SM 



2ir 



SMi 
12 I 



the quantity Tf^ 1 is given by ]16 



^SM 
12 



v t + - v c v t \z c + -z c - -z c + 



G\ m\ M Bd B Bd p Bd 
8tt 

n 3— A 2 \ 8 2 2 1 , 



(17) 



3 3 c 3 c 

where Vi = V^, V* d and z c = m 2 c jm\. Varying f Bd , \V u b\, and \V c b\ in the range quoted in 
Table |, and 7 inside the range given by (54.8 ± 6.2)° fl7| , we get the SM value to be 

-1.54 x 10~ 3 < A u 

whereas the current world average is M 



SM < -0.16 x 10~ 3 , 



A^ p « (0.2 ± 1.4) x 10- 2 . 

In the presence of SUSY, the phases of Mi 2 and Ti 2 may be no longer the same, and 
potentially a larger dilepton charge asymmetry may be possible. In particular, M 12 could 
be affected strongly by SUSY particles, whereas T 12 is not, i.e. Af£ ULL = Mf 2 M + Mf 2 USY 
whereas T™ 1 " 1 " « l~j 2 . In this case, the dilepton charge asymmetry could be approximated 

as 



A, 



Im 



1 12 



Mf 2 M + MS, USY 



12 



(18) 



The possible ranges of An in a class of general SUSY models were studied in Ref. [18 
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B. Effective Hamiltonian for AB = 1 processes 



The effective Hamiltonian relevant to AB = 1 processes involves four quark operators 
and b — > and b — > dg penguin operators. Since we are not going to discuss AB = 
1 nonleptonic decays due to theoretical uncertainties related with factorization, we shall 
consider the inclusive radiative decay B — ► X d ^ only The relevant effective Hamiltonian for 
this process is given by [19| 



H e s(b -> dr/(+g)) = -j=-V t * d V tb \^ Ci{fi b )O ic {^ b 

V ^ 1=1,2.7,8 



+ 



AC 

-&V: d V ub J2 Pi»(A*6) - O fc (A*6)] (19) 

4=1,2 



with 



n 



Oic = ^l7 m cl cr7 M 6x, O iu = d L ^u L u L j,J) L , 
2c = di/fc L clj^l, 2u = di/fu L ulI^l, 

' m b d L a^F„ v b R , 8a = ^m b d L a^t a G^b R - 



lQn 2 



16vr 2 



(20) 



Here the renormalization scale fi b is of the order of m b , and we have used the unitarity of 
the CKM matrix elements 



v; d v cb = -(v: d v ub + v td v t 



u v tb), 

which should be valid even in the presence of SUSY flavor violations. 

In the SM, all the three up-type quarks contribute to this decay, since all the relevant 
CKM factors are of the same order of magnitude. The strong phases are provided by the 
imaginary parts of one loop diagrams at the order 0(a s ) by the usual unitarity argument. 
Varying f Bd , \V ub \, and \V cb \ in the range quoted in Table |, and 7 between (54.8 ± 6.2)° 



|17| , we get the branching ratio for this decay in the SM to be 8.9 x 10~ 6 — 1.1 x 10 -5 . The 

direct CP asymmetry in the SM is about —15% 10% [[HJ. We have updated the previous 

predictions by Ali et al. |19 using the present values of CKM parameters. 

The CP averaged branching ratio for B — ► X d ^ in the leading log approximation is given 



by [19, 20, 21 



B(B -> X dl ) 
B(B — > X c eu) 



v td v tb 



cb 



Qa 



|C 7 (m fe )p 



(21) 



where f(z) = 1 — 82; + 8^ 3 — z A — 12z 2 In z is the phase space factor for the b — > c semileptonic 
decays and a -1 = 137.036. Neglecting the RG running between the heavy SUSY particles 
and the top quark mass scale, we get the following relations : 



C 7 (m b ) 
C 8 (m b ) 



-0.31 + 0.67 CI 



mew 

'8 



(m w ) + 0.09 C, 
(m w ). 



new 
8 



-0.15 + 0.70 CI 

The new physics contributions to C 2 are negligible so that we use C 2 {m b ) 



Ci M (m b ) 



(22) 
1.11. 



7 



u 8g 



In general SUSY models considered in the present work, the Wilson coefficients G> 

1 



are given by M. 15 



7 n 7 and 



C™ SY (m w ) 



8nQ b a s 



3V2G F mW t * d V tb _ 
27TQ;, 



V2G F m%*V t 



(5 



(5f 3 ) LL M 4 (x 
3 



( S 13) LR 



m\ x 



m b 



ABAx) 



(23) 



11) LL 



+(5 d 



13) LR 



td v tb L 

rh\/x 



2 M 3 (x) - -M 4 (x) 



- (AB l (x)-9x- 1 B 2 (x)) 



(24) 



Here we have ignored the RG running between the squark and the gluino mass and the 
mw scale. Note that the {Sf 3 )LR contribution is enhanced by mg/m b compared to the 
contributions from the SM and the LL insertion due to the chirality flip from the internal 
gluino propagator in the loop. Explicit expressions for the loop functions B^s and Mj's can 
be found in Ref. 



is, m. 



In order to generate a nonvanishing direct CP asymmetry, one needs at least two indepen- 
dent amplitudes with different strong (CP-even) and weak (CP-odd) phases. In B — > Xdj, 
strong phases are provided by quark and gluon loop diagrams, whereas weak phases are 
provided by the KM angles (a, (3, 7) and (Sf 3 )AB- The resulting direct CP asymmetry in 
B — > X d j can be written as Jl9|, |20| 



4p^(in 



[10.57 Im (C 2 C 7 *) - 9.40 Im ((1 + e d )C 2 C* 7 



|CV. 

9.51 Im (C 8 C 7 *) + 0.12 Im ((1 + e d )C 2 C* 8 



(25) 



where 



V td V tb ~ (1-p + iv) 



in the Wolfenstein parametrization for the CKM matrix elements. 

A remark is in order for the above CP asymmetry in B — > Xd'j- Unlike the B 



Xg'j case 



for which the |CV 7 | is constrained by the observed B — > X s ^/ branching ratio, the B — > Xd'j 
decay has not been observed yet, and its branching ratio can be vanishingly small even in 
the presence of new physics. In that case, |C- 1 n °" ' , '"" ,, " ;, — f " v " r 



77 1 



so that the denominator of A 



CP 



becomes zero and the CP asymmetry blows up. This could be partly cured by replacing the 
denominator |C 77 | 2 by i^NLo(^) defined in Ref. p0| : 



0.11 \C 2 \ 2 + 68.13 \C 7 \ 2 + 0.53 |C 8 | 2 - 16.55 Re(C 2 C; 



- 0.01 Re(C 2 C;) + 8.85 Re(C 7 C 8 *) + 3.86 Re(C{ 1) C 7 * 



(26) 



for the photon energy cutoff factor 5 
tion to C 7l {m b ) pOf : 



0.3. Here C 7 is the next-to-leading order contribu- 



°77 



0.48 - 2.29 C^{m w ) - 0.12 C£ cw (m w ). 
This prescription will render the denominator of (|25|) to be finite. 



(27) 
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III. NUMERICAL ANALYSIS 



In the numerical analysis, we impose the following quantities as constraints : 

• Anifl^ = (0.472 ± 0.017) ps _1 : This constrains the modulus of Mi2(B°) through the 
following relation Am Bd = 2 \M 12 (B°)\ §. 

• = (0.79 ± 0.10) : This constrains the phase 20 of M 12 (B ) by = sin 20, 
where 2f3 is the argument of Mi2(-B°) 0. 

• BriB — > AVy) < 1 x 10~ 5 : At present, there are limits only on the exclusive decays: 
BiB — > /ry) < 2.3 x 10~ 6 . We assume a modest upper bound on the branching ratio 
for the inclusive radiative decay as BriB — > Xffi) <lx 10~ 5 . This is also well below 
the experimental uncertainty in the B — > branching ratio. This puts a strong 
constraint on both LL and LR insertions as we shall see. Especially the LR insertion 
is more strongly constrained, since its contribution is enhanced by m g /mi, due to the 
chirality flip from the gluino in the loop compared to other contributions including 
the SM one. This is a new ingredient compared to the work in Ref. f2|. 

• A u xp = (0.2 ± 1.4)%: This is related to the B° — B° mixing through the relation (16). 
Although we do not use this constraint to restrict the allowed parameter space, we 
indicate the parameter space where the resulting An falls out of the lcr range. It turns 
out that both LL and LR mixing scenarios are already strongly constrained by the 
B — > Xj'j branching ratio rather than by An, as can be seen in Figs. 2 (a) and (b). 

We impose these constraints at 68 % C.L. (1<t) as we vary the KM angle 7 between and 2ir. 
In all cases, we set the common squark mass fh = 500 GeV and x = 1 (trig = fh). Finally 
for the mass insertion parameters {^ 3 )ab, we consider two cases. In the first case (the LL 
case), only (5f 3 )iL is nonvanishing among the mass insertion parameters, and in the second 
(the LR case), only (5f 3 )LR- It would be straightforward to consider other possibilities such 
as (^i 3 )ll = (8f 3 )fiR etc. in a similar way. 

The parameter space searching is done in the following way. We vary 7 from to 2tt, and 
(^13) ab inside the bound depicted in Ref. @. For a given set of values of 7 and (6f 3 )AB, we 
search for those fg d and y/ p 2 + rf = \V u b\/XV c b (with A = \V US \) that satisfy 1-er constraints 
on AMb and sin 20 . The search region is the 1-a range in Table |I[ If no such pair exists, 
the (7, (#13) ab) point is excluded from the plots. Points that are not excluded are plotted in 
Fig. |l|. Using these 7, (5f 3 )AB, fB d , and \J p 2 + rf found above, we plot BriB — > X^j) and 
A 1 ^ 1 . Uncertainties in .Bi_ 3 (m&), which are actually used in our analysis, are only several 
per cent, while that in fs d is 15%. Moreover the matrix elements flllD are proportional to 
f B Bi(mb), so we do not take into account the uncertainties in the bag parameters. 

In Figs. (a) and (b), we show the allowed parameter space in the (Re(5f 3 )AB, Ihi^^ab) 
plane [(a) LL insertion and (b) LR insertion, respectively] for different values of the KM 
angle 7 with different color codes: dark (red) for 0° < 7 < 90°, light gray (green) for 90° < 
7 < 180°, very dark (blue) for 180° < 7 < 270° and gray (magenta) for 270° < 7 < 360°. 
The region leading to a too large branching ratio for Bd — > AVy is covered by slanted lines. 
And the region where An falls out of the data within la range is already excluded by the 
B — > Xd'j branching ratio constraint. For both the LL and LR mixing cases, our results are 
the same as those in Ref. 0, if we impose only the Ams d and sin 2/3 constraints. By adding a 
constraint from B — > AVy (and An), the allowed parameter space is further reduced, and the 
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-0.3 




-0.08 



_i_ 



_i_ 



_i_ 



-0.15 



0.15 



Re (5? 3 )ll 
(a) LL mixing only 



-0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 

Re (5? 3 )lr 

(b) LR mixing only 



FIG. 1: The allowed ranges in (a) the LL and (b) the LR insertion cases for the parameters 
(Re(5f 3 )AB-, I m (^i3)A_B) for different values of the KM angle 7 with different color codes: dark 
(red) for 0° < 7 < 90°, light gray (green) for 90° < 7 < 180°, very dark (blue) for 180° < 7 < 270° 
and gray (magenta) for 270° < 7 < 360°. The region leading to a too large branching ratio for 
Bfi — > Xjpf is colored lightly and covered by parallel lines. 



~i 1 r 

b->dy allowed 




50 100 150 200 250 300 350 

y n 



(a) LL mixing only 



(b) LR mixing only 



FIG. 2: The possible ranges of the dilepton charge asymmetry in (a) the LL and (b) the LR cases 
as functions of the KM angles 7. The black rectangle around 7 ~ 55° is the SM prediction. Those 
parameters which lead to B(B — > X^j) > 1 x 1CP 5 are denoted by the the gray (magenta) region, 
and those for B(B — > X^i) < 1 x 1CP 5 by the dark (blue) region. The la range for the world 
average of A^ p = (0.2 ± 1.4)% is shown to lie between the short dashed lines. 
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FIG. 3: The possible ranges of (a) B(Bd — > Xd"i) and (b) A^ rf7 as functions of the KM angle 7 in 
the LL insertion case. The black rectangle around 7 ~ 55° is the SM prediction. Those parameters 
which lead to B(B — ► X^p/) > 1 x 10~ 5 are represented by the gray (magenta) region, and those 
for B(B —> Xdj) < 1 x 10~ 5 by the dark (blue) region. 

effect is even larger for the LR mixing case. For the LL mixing [ Fig. 1 (a) ], the B — > X<f] 
does play some role, and the An gives a moderate constraint. The KM angle 7 should be 
in the range between ~ —60° and ~ +60°, and An can have the opposite sign compared to 
the SM prediction, even if the KM angle is the same as its SM value 7 ~ 55°. For the LR 
mixing [ Fig. 1(b)], the B(Bd — > Xd'j) puts an even stronger constraint on the LR insertion, 
whereas the An does not play any role. In particular, the KM angle 7 can not be too much 
different from the SM value in the LR mixing case, once the B(Bd — > X^y) constraint is 
included. Only 30° < 7 < 80° is compatible with all the data from the B system, even if we 
do not consider the ex constraint. The resulting parameter space is significantly reduced 
compared to the result obtained in Ref. 0. The limit on the LR insertion parameter will 
become even stronger as the experimental limit on Bd — > Xd'j will be improved in the future. 

In Fig. |^, we show the predictions for An as a function of the KM angle 7: (a) LL insertion 
and (b) LR insertion only. For the LL insertion case [Fig. 2 (a)], one can expect a large 
deviation from the SM prediction for An for a wide range of the KM angle 7 (~ —60° < 7 < 
+60°), even after we impose the B — > Xd'j branching ratio which is more constraining than 
the An (the short dashed lines indicate the la range for A u xp ). Also even if the KM angle 
7 is close to the SM value (7 ~ 55°), the dilepton charge asymmetry An can be different 
from the SM prediction by a significant amount due to the SUSY contributions from (5^ 3 )ll. 
On the other hand, for the LR insertion case [Fig. 2 (b)], the B — > Xdj constraint rules 
out essentially almost all the parameter space region, and the resulting An is essentially the 
same as for the SM case. 

In Figs. (a) and (b), we show the branching ratio of Bd — > Xd'J and the direct CP 
asymmetry therein, respectively, as functions of the KM angles 7 for the LL insertion only. 
The SM predictions 

B{B d -> Xdi) = (0.9 - 1.1) x 1(T 5 , A b ^ dl = (-15 ~ -10)% 
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FIG. 4: The possible ranges of (a) B(Bd — > and (b) A^^ 7 as functions of the KM angle 

7 in the LR insertion case . The black rectangle around 7 ~ 55° is the SM prediction. Those 
parameters which lead to B(B — > ^7) > 1 x 10 -5 are represented by the gray (magenta) region, 
and those for B(B — > ^7) < 1 x 10~ 5 by the dark (blue) region. 



are indicated by the black boxes. In this case, the KM angle 7 is constrained in the range 
~ —60° and ~ +60°. The direct CP asymmetry is predicted to be between ~ —15% and 
~ +20%. In the LL mixing case, the SM gives the dominant contribution to Bd — > Xd'y, 
but the KM angle can be different from the SM case, because SUSY contributions to the 
B° — B° mixing can be significant and the preferred value of 7 can change from the SM KM 
fitting. This is the same in the rare kaon decays and the results obtained in Ref. [15| apply 
without modifications. If the KM angle 7 is substantially different from the SM value (say, 
7 = 0), we could anticipate large deviations in the Bd — > Xdj branching ratio and the direct 
CP violation thereof. 

In Figs. |4] (a) and (b), we show the branching ratio of Bd — > Xd'y and the direct CP 
asymmetry therein, respectively, as functions of the KM angles 7 for the LR insertion only. 
As before, the black boxes represent the SM predictions for B{Bd —* Xdj) and the direct CP 
asymmetry therein. In the LR insertion case, there could be substantial deviations in both 
the branching ratio and the CP asymmetry from the SM predictions, even if the Am# and 
sin 2/5 is the same as the SM predictions as well as the data. For the LL insertion, such a 
large deviation is possible, since the KM angle 7 can be substantially different from the SM 
value. On the other hand, for the LR mixing, the large deviation comes from the complex 
(^13) lr even if the KM angle is set to the same value as in the SM. The size of {S( 3 )lr, is 
too small to affect the B° - W mixing, but is still large enough too affect B — > Xdj- Our 
model independent study indicates that the current data on the Amg, sin 2/? and An do 
still allow a possibility for large deviations in B — > X^7, both in the branching ratio and the 
direct CP asymmetry thereof. The latter variables are indispensable to test completely the 
KM paradigm for CP violation and get ideas on possible new physics with new flavor/CP 
violation in b — > d transition. 
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IV. CONCLUSIONS 



In this work, we considered the gluino-mediated SUSY contributions to B° — B° mixing, 
B — > J/ipK s and B — > in the mass insertion approximation. We find that the (I/L) 
mixing parameter can be as large as | (^13) ii| ^ 2 x 10 _1 , but the (LR) mixing is strongly 
constrained by the B — > X^7 branching ratio: | (^13) z,/? I ^ 10~ 2 . The implications for 
the direct CP asymmetry in B —>■ X^ are also discussed, where substantial deviations 
from the SM predictions are possible both in the LL and LR insertion cases for different 
reasons. For the LL insertion case, the SUSY contribution to B — > is not so significant, 
but is still constrained by the current upper limit on B — > Xd'j- (If the upper limit were 
B(B -> X d >y) < 5 x 10" 5 , then the allowed region for the KM angle 7 is the whole range from 
to 2tt, as can be seen from Fig. 3 (a). In this case, the An will provide a more important 
constraint for the LL insertion.) Also the global KM fitting can change because SUSY can 
affect B° — B° mixing by a significant manner. By the same reason, there is still ample room 
for large deviations in the An for the LL insertion case. On the other hand, for the LR 
insertion case, the SUSY contribution to B — > X^ is enhanced by the factor m g /mb and 
the size of (Sf 3 )iR is strongly constrained. The resulting effect is that the KM angle cannot 
differ too much from the SM case. Still large deviations in the branching ratio for B —>■ Xa'y 
and direct CP violation thereof is possible due to large SUSY loop effects on B — > X d j. 
Thus it turns out that all the observables, An, the branching ratio of B — > X^ and the 
direct CP violation thereof are very important, since they could provide informations on 
new flavor and CP violation from (b~fz)LL,LR (or any other new physics scenarios with new 
flavor/CP violations). Also they are indispensable in order that we can ultimately test the 
KM paradigm for CP vioaltion in the SM. 

Note Added 

While this work was being finished, we recieved a preprint , in which similar processes 
(the exclusive B — > /ry and various asymmetries thereof, and An) in certain class of SUSY 
models are considered. 



Acknowledgments 

This work is supported in part by BK21 Haeksim program of the Ministry of Education 
(MOE), by the Korea Science and Engineering Foundation (KOSEF) through Center for 
High Energy Physics (CHEP) at Kyungpook National University, and by DFG-KOSEF 
Collaboration program (2000) under the contract 20005-111-02-2 (KOSEF) and 446 KOR- 
113/137/0-1 (DFG). 



[1] B. Aubert et al. [BABAR Collaboration], fcrXiv:hep-ex/0203007| T. Higuchi [Belle Collabo- 
ration], |arXiv:hep-ex/020502q . 



[2] D. Becirevic et al, |arXiv:hep-ph/0U2303| . 

[3] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 
[4] S. Laplace, Z. Ligeti, Y. Nir and G. Perez, Phys. Rev. D 65, 094040 (2002) |arXiv:hep- 
ph/020201CH . 



13 



[5] A. Ali, G. F. Giudice and T. Mannel, Z. Phys. C 67, 417 (1995) [|arXiv:hep-ph/9408213|l ; 
Y. G. Kim, P. Ko and J. S. Lee, Nucl. Phys. B 544, 64 (1999) larXiv:hep-ph/9810336| ; A. Ali, 
E. Lunghi, C. Greub and G. Hiller, |arXiv:hep-ph/011230q . 

[6] F. Gabbiani and A. Masiero, Nucl. Phys. B 322 (1989) 235 ; J. S. Hagelin, S. Kelley and 
T. Tanaka, Nucl. Phys. B 415 (1994) 293 ; E. Gabrielli, A. Masiero and L. Silvestrini, Phys. 
Lett. B 374 (1996) 80 jarXiv:hep-ph /9509379| ; F. Gabbiani, E. Gabrielli, A. Masiero and 
L. Silvestrini, Nucl. Phys. B 477 (1996) 321 [ |arXiv:hep-ph/96043871 ; 

[7] A. Ali and E. Lunghi, Eur. Phys. J. C 21, 683 (2001) [|arXiv:hep-ph/010520qi ; A. Ali, Contri- 
bution to International Conference on Flavor Physics (ICFP 2001), Zhang- Jia-Jie City, Hunan, 
China, 31 May - 6 Jun 2001. |arXiv:hep-ph/020112q . 

[8] L. Everett, G. L. Kane, S. Rigolin, L. T. Wang and T. T. Wang, JHEP 0201, 022 (2002) 
[|arXiv:hep-ph/0112T26| . 

[9] S. Baek and P. Ko, Phys. Rev. Lett. 83, 488 (1999) [|arXiv:hep-ph/9812229|| ; A. Ali and 

D. London, Eur. Phys. J. C 9, 687 (1999) ||arXiv:hep-ph/9903535|| ; A. Bartl, T. Gajdosik, 

E. Lunghi, A. Masiero, W. Porod, H. Stremnitzer and O. Vives, Phys. Rev. D 64, 076009 
(2001) [|arXiv:hep-ph/0103324|l ; A. J. Buras and R. Fleischer, Phys. Rev. D 64, 115010 (2001) 
[|arXiv:hep-ph /0104"238| . 

[10] M. Convery [BABAR Collaboration], talk presented at the meeting of the APS Division of 

Particles and Fields DPF-2002, Williamsburg, Vriginia, 24-28 May, 2002. 
[11] A. J. Buras, in 'Probing the Standard Model of Particle Interactions', F.David and R. Gupta, 



eds., 1998, Elsevier Science B.V. |arXiv:hep-ph/9806471 



[12] A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B 347, 491 (1990). 

[13] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) |arXiv:hep- 

| ph/9512380H - 

[14] D. Becirevic, V. Gimenez, G. Martinelli, M. Papinutto and J. Reyes, JHEP 0204, 025 (2002) 
|arXiv:hep-lat/0110'091 1. 

[15] S. Baek, J. H. Jang, P. Ko and J.H. Park, Nucl. Phys. B 609, 442 (2001) |arXiv:hep 



ph/0105028;1 



[16] A. J. Buras and R. Fleischer, Adv. Ser. Direct. High Energy Phys. 15, 65 (1998) | arXiv:hep- 



ph/970437l . 

[17] M. Ciuchini et al, JHEP 0107, 013 (2001) |arXiv:hep-ph/0012308 |. 
[18] L. Randall and S. Su, Nucl. Phys. B 540, 37 (1999) [|arXiv:hep-ph/9807377| ; 
[19] A. Ali, H. Asatrian and C. Greub, Phys. Lett. B 429, 87 (1998) [bxXiv:hep-ph/9803314 ]. 
[20] A. L. Kagan and M. Neubert, Phys. Rev. D 58, 094012 (1998) [ |arXiv:hep-ph/98033681 . 
[21] A. G. Akeroyd, Y. Y. Keum and S. Recksiegel, Phys. Lett. B 507, 252 (2001) [|arXiv:hep- 
ph/010300^]. 



[22] A. J. Buras, G. Colangelo, G. Isidori, A. Romanino and L. Silvestrini, Nucl. Phys. B 566, 3 

(2000) jarXiv:hep-ph/99083711 ; T. Besmer, C. Greub and T. Hurth, Nucl. Phys. B 609, 359 

(2001) [|arXiv:hep-ph/0105292|| ; 

[23] A. Ali and E. Lunghi, |arXiv:hep-ph/0206242| . 



14 



